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The data for 9.3 miUion T{2S) and 20.9 milhon T(15') taken with the CLEO III detector has been 
used to study the radiative population of states identified by their decay into twenty six different 
exclusive hadronic final states. In the T(25') decays an enhancement is observed at a ~ 5(t level 
at a mass of 9974.6 ± 2.3(stat) ± 2.1(syst) MeV. It is attributed to rib{2S), and corresponds to the 
T{2S) hyperfine splitting of 48.7± 2.3(stat) ± 2.1(syst) MeV. In the T(15') decays, the identification 
of ryi,(15) is confirmed at a Stj level with M(j7(,(15)) in agreement with its known value. 



PACS numbers: 14.40.Pq, 12.38.Qk, 13.25.Gv 

The spin-dependent interaction between constituents 
of a composite system is of general interest, from hy- 
drogen and positronium to mesons and baryons. The 
hyperfine interaction between quarks is one of the most 
important components of the spin-dependent QCD inter- 
action. It has been studied for charm quarks by measure- 
ments of hyperfine sphttings between spin-triplet and 
spin-singlet states, M{n^L) — M{n}L), of IS, 2S, and 
IP states of charmonium [1]. For the bottom quarks, 
measurements have been reported recently for the IS hy- 
perfine splittings [2] and for the IP and 2P hyperfine 
splittings [3] of bottomonium states. While ?7f,(lS'), the 
bottomonium ground state, has been identified, the radi- 
ally excited r]i,{2S) has not been identified so far, and 
the hyperfine splitting of the bottomonium 2S states, 
AMhf{2S) = T{2S) - r]b(2S), is not known. In this 
letter, we report on the observation of r]ii{2S) in its for- 
mation in the radiative decay of T{2S) and its exclusive 
decays into twenty-six different final states containing 
charged light-quark hadrons, pion, kaons, and antipro- 
tons. A similar analysis of T(IS') decays is made, and it 
confirms the identification of ?7f,(lS'). 

An early attempt to identify r/b(lS', 25) in the inclusive 
allowed Ml radiative decays, T(15,2S') -frib{lS,2S) 
by detecting the low energy (< 100 MeV) transition 
photons was unsuccessful [4] . Successful identification of 
77f,(lS') by BaBar [2], and its confirmation by CLEO [2], 
was only made possible by detecting the ~ 920 MeV 
and ~ 610 MeV transitions photons in the "forbidden" 
Ml decays T{3S,2S) IVbi^S). These Ml transitions 
have zero overlap between the initial and final states 
in the lowest order. They become finite only because 
of relativistic and higher-order effects, and theoretical 
predictions for them are notoriously difficult and unre- 
liable. For example, a recent calculation in the frame- 
work of non-relativistic effective field theory of QCD (pN- 
RQCD) [5] predicted a branching fraction for the decay 
T(2S') -^' 7776(15) more than two orders of magnitude 
larger than that measured by BaBar [2]. In contrast, 
"allowed" Ml transitions, T{nS) ^rjb{nS) are rela- 
tively simple, with the wave function overlaps between 



the initial and final states being essentially unity. They 
are therefore very attractive from the theoretical point of 
view. Since the inclusive radiative transitions are essen- 
tially impossible to measure, the only hope is to identify 
r]b{nS) in radiative decays of T{nS) by "tagging" r]b{nS) 
by their exclusive hadronic decays. 

In this paper, we report on the study of the reaction 

T(25) -> ^Vb{2S), 77,(25) ^ X, 

(X=4,6,8,10 7r±,if±,p/p) 

We also report on an identical study of T(15) 
7?7f,(15), 77;,(15) X, which provides a useful check of 
our analysis procedure. We use data obtained with the 
CLEO III detector at the Cornell Electron Storage Ring, 
CESR. The data consist of (9.32 ±0.19) x 10^ T(25) and 
(20.82 ± 0.41) X 10^ T(15). To develop event selection 
criteria we use Monte Carlo (MC) simulations and data 
consisting of (5.88 ± 0.12) x 10^ T(35). 

The CLEO III detector, which has been described be- 
fore [6], consists of a Csl electromagnetic calorimeter, an 
inner silicon vertex detector, a central drift chamber, and 
a ring-imaging Cherenkov (RICH) detector, all inside a 
superconducting solenoid magnet with a 1.5 T magnetic 
field. The detector has a total acceptance of 93% of Att 
for charged and neutral particles. The photon energy res- 
olution in the central (83% of 47r) part of the calorimeter 
is about 2% at = 1 GeV and about 5% at 100 MeV. 
The charged particle momentum resolution is about 0.6% 
at 1 GeV/c. 

We select events that have 4, 6, 8, or 10 charged par- 
ticle tracks with zero net charge and at least one pho- 
ton candidate. The charged tracks are required to be 
well-measured and consistent with coming from the in- 
teraction point. The photon candidates are calorimeter 
showers which lie within the "good barrel" or "good end- 
cap" regions, |cos0| < 0.81 and |cos(?| ~ 0.85 — 0.93, 
respectively, where is the polar angle with respect to 
the incoming positron direction. They are required to 
contain at least 10 MeV of energy, to not contain any 
of the few known noisy calorimeter cells, and to have 
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a transverse energy distribution consistent with an elec- 
tromagnetic shower. Analysis of simulated events shows 
that the largest background to the low energy transition 
photons comes from calorimeter showers due to the in- 
teraction of the final state charged hadrons with detector 
elements. To reduce this background, we make the isola- 
tion requirement that a photon candidate must be sepa- 
rated from the nearest charged track by 50 cm. To ensure 
that there is no contribution from photons from tt^ — > 77 
decays, we reject photon candidates which make a two- 
photon invariant mass of M(77) = M(7r'^)±25 MeV with 
any other photon candidate in the event. 

Charged tracks are identified as tt^, K^, and p/p us- 
ing dE/dx, the energy loss in the drift chamber, and 
information from the RICH subdetector. To utilize 
dE/dx information, for each particle hypothesis, X = 

TT, K, p or p, we calculate Xx^^^ = [(d-E/dx)measured - 

(di^/fix) predicted] /cx, for hypothcsis X, and ax is the 
standard deviation of the measured dE/dx for hypothe- 
sis X. For higher momentum tracks with | cos 6*1 < 0.8, 
we use the combined log-likclihood variable 
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where tRICH 

pothesis obtained from measurement in the RICH sub- 
detector. 

For low momentum tracks with p < 0.6 GeV for 
TT* and K^, and p < 1.5 GeV for p/p, we only use 
dE/dx information, and require the measured dE/dx to 
be within 3cr of the expected energy loss for the parti- 
cle hypothesis. We also require |xit'^^'''^| < Ix^'^'^^l for 
charged pions, Ix^^''^! < |Xit^^''^I for charged kaons, 
and Ixp^^'^^l < \x'^K^'^^\ for p/p. For higher momen- 
tum tracks, we require ^Lt^^k < for charged pions, 
ACk,7i < for charged kaons, and ACp,K < for p/p. 
We remove contamination of electrons by rejecting events 
with tracks which have a ratio between the energy de- 
posited in the calorimeter and the momentum measured 
in the drift chamber, E/p — 0.9 — 1.1. We reconstruct 
Kg — >■ 7r+7r~ by requiring these decays to have a common 
vertex that is displaced > 3a from the interaction point. 

We reconstruct the f]b{nS) candidate in the following 
26 decay modes: 2(7r+7r~), 3(7r+7r" 
K+K^TT+TT-, K+K-2{tt+tt-) 



), 4(7r+7r-), 5(7r+7r-) 

K+K-'i{'K+TT^) 



K+K-4:{tt+tt-), 

pp2(7r+7r^), pp3{'K' 
ppK+K-2{TT+n 



2{K+K-), 
2{K+ K-)?,{tt^ 



ppK+K-2{TT+n~), ppK+K-^-K+TT-), K^K'^ir^ 
if^i^±7r=F7r+7r-, K'^K^Tr^2(Tr+Tr-), K K^ir"^ 3 (tt+tt-) 
2KIti+ti-, 2i^g2(7r+7r-), 2i^g3(7r+7r-). 

To select events with well-measured hadrons, we fit the 
reconstructed hadrons to a common vertex, and require 
that the reduced x^ of the fit is /d.o.f. < 4. 

To reconstruct the full event including both the 
hadrons and the transition photon, we perform a 4C 
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FIG. 1. Invariant mass distributions in terms of AA/ = 
M{r{2S, IS)) - M(hadrons) for (top) T(2S') and (bottom) 
T(15') data with final event selections. 



kinematic fit constraining the combination of a photon 
candidate and the hadronic final state to have the center- 
of-mass four-momentum of zero (except for a small con- 
tribution due to the finite beam crossing angle) . We per- 
form this fit for each signal photon candidate in the event 
and pick the fit with the lowest x^- Henceforth, we use 
the constrained hadronic mass which has the much bet- 
ter mass resolution (a ~ 5 MeV) of the photon. To 
reject T{nS) hadrons events combined with a fake 
photon, we use the reduced x^ of the 4C fit. We require 
X^/d.o.f. < 4. 

As noted by BaBar, because there is no preferred di- 
rection in the decay of the spin-zero ry^, there is weak 
correlation between the signal photon momentum in the 
center-of-mass frame with the thrust axis calculated for 
the hadrons from the decay of 77b. In contrast, the same 
correlation is strong for the background events. There- 
fore, the signal-to-background ratio varies with the angle 
6t between the photon and the thrust vector, and a cut 
on |cos0t| is very useful in rejecting background. This 
was confirmed by both BaBar and CLEO in their iden- 
tification of rib{lS), and is found to be also true in MC 
simulations in the present case. A cut to accept events 
with I cos^tI < 0.5 is found to be optimum. 

The efficiencies for individual decay modes were de- 
termined by MC simulations. They range from 8.2% to 
0.5%, depending on event multiplicity. 

The invariant mass distributions for T(IS') and T(2S') 
data with final event selections are shown in Fig. 1 in 
terms of AM = M{T{nS)) - M(hadrons). At the 
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FIG. 2. Distributions of AM = M{T{2S, IS)) - A'/(hadrons): (top row) T(25') data, (bottom row) T(1S) data are shown with 
both (left column) linear and (right column) logarithmic scales. The best fit curves are shown as the thick sohd curve and the 
individual components are shown as the dashed and thin solid curves. 



TABLE I. Results of fits to the AA'I{2S, IS) data distributions, as described in the text. Only statistical uncertainties are 
given. The product branching fraction is 81x82= Bi(T{nS) — >■ 7?7(,(nS)) x J2i=i ^'^iiVbinS) hi). The first errors are 
statistical, and the second errors are systematic, as detailed in Table II. 





N 


AMhf (MeV) 


PI (MeV) 


x'/d.o.f. 


signif. (a) 


B1XB2X 10" 


rib{2S) 
Vb{lS) 


llA+ti 


48.7 ±2.3 ±2.1 
67.1 ±3.4 ±2.3 


9974.6 ± 2.3 ± 2.1 
9393.2 ±3.4 ±2.3 


91.8/103 
114.6/107 


4.9 
3.1 


46.2+^^;^ ± 10.6 
SO.ltr^f ± 7.5 



smallest values of AM, the distributions show the large, 
steeply falling contributions due to T{nS) decays. En- 
hancements are seen at AM « 50 MeV in the T{2S) data 
and at AM w 70 MeV in the T(IS') data. These appear 
to be well separated from the T contributions. In the 
T(2S') data, the contributions from T{2S) JXbji^P), 
Xbj{^P) — >■ hadrons are also seen in the range AM w 
100 - 200 MeV. The distributions in Fig. 1 have been 
fitted as follows. 

The non-peaking backgrounds due to continuum 
and misidentifications in the T{2S) data is essen- 
tially zero both below and above the Xbj{iP) peaks, 
and in the T{1S) data it is nearly constant (^ 
0.4 count/2.5 MeV bin) in the region AM = 100 - 
300 MeV. 

The choice of the fit function for the large, rapidly 
falling T(2S') and T(IS') contributions at small AM is 
important. It is not possible to obtain the total shapes 
of the T(2iS') and T(IS') contributions from MC simula- 
tions. We have made a large number of MC simulations, 
and find that the predicted shapes of these contributions 
differ substantially between different decays and multi- 
plicities, particularly in the tail regions. Further, they 



can not be added to produce the composite shape be- 
cause the relative proportions of the individual contribu- 
tions are not known. Hence, an empirical approach to fit 
it was adopted, and fits with different fit functions (ex- 
ponentials of the form exp(aa; + bx'^ + cx^ + ■■")) were 
tried. The best fits were consistently obtained with a 
single exponential. Single exponentials were also found 
to best fit the data for T(3S'). In the left panels of Fig. 2 
we show the fits in linear plots. In the right panels we 
show the same fits in log plots to illustrate that the single 
exponentials fit the T{nS) contributions very well, and 
the enhancements at ~ 70 MeV and ~ 50 MeV in the 
T(IS') and T{2S) data received very little contribution 
from the "tails" of T(15,25'). 

We fit the peaks in the AM distributions with Breit- 
Wigner shapes convolved with the known Gaussian ex- 
perimental resolution functions which have widths which 
vary from a = 4.2 MeV at AM = 50 MeV to ct = 
6.4 MeV at AM = 165 MeV. The Breit-Wigner width of 
the enhancement in T(25) at AM w 50 MeV is assumed 
to be 5 MeV. The Breit-Wigner width of the enhance- 
ment in T(IS') at AM w 70 MeV, attributed to r/b(lS'), 
is assumed to be 10 MeV. The Xbj{^P) peaks are fitted 
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with the Gaussian rcsohition widths. The masses of the 
Xbj{^P) peaks are found to be in agreement with their 
known masses within 1.1 ± 0.8 MeV on average. 

The fit results are listed in Table I. 

The significance values in Table I are determined as 
a = 

max); where Lmax is the maximum 
likelihood returned by the fits including the enhance- 
ments at ~ 70 MeV and ~ 50 MeV, and Lq is the hkeli- 
hood returned by fits without these enhancements. 

The fitted value of the ^ 70 McV enhancement in 
the T(15) data is AM = 67.1 ± 3.4(stat) MeV, and 
the observation has a significance of 3. la. It is nat- 
urally identified as being due to 77f,(lS'), and leads to 
AMhf (15)65 = 67.1 ± 3.4(stat) MeV, in good agreement 
with the PDG average of 69.3 ± 2.8 MeV [1]. The fitted 
value of the ^ 50 MeV enhancement in the T(25') data 
is 48.7 ± 2.3(stat) MeV, and the observation has a sig- 
nificance of 4.9(7. We can not find any explanation for 
this 4.9(7 enhancement except to attribute it to rib{2S). 
Henceforth, we refer to it as such. Thus, we determine 
AMhf (2S')f,5 = 48.7±2.3(stat) MeV. This constitutes the 
first determination of hyperfine splitting in the bottomo- 
nium 25" radial excitation. 

To confirm the likelihood determination of the signif- 
icance of the rj}j{2S) enhancement, we have made a MC 
determination (taking systematic uncertainties into ac- 
count) of the probability that the observed enhancement 
can arise due to a statistical fluctuation anywhere in the 
range AM = 35 - 70 MeV. In 10^ trials, 4565 such fluc- 
tuations were found. This corresponds to a significance 
of 4.6cr for the observed enhancement. 

We calculate product branching fractions correspond- 
ing to the observed counts A'^i(obs) in individual hadronic 
decay modes hi using MC-determined efficiencies as 



TABLE II. Summary of systematic uncertainties and their 
sums in quadrature. 



Bi[T(n5)^7%(n5),7X6j(lJ')] 

X S2j[(??6,X6j) hi] = 



N^{ohs) 



ti X N{T{nS)) ' 



(1) 



where N {T (nS)) refers to the number of T(nS) in the 
data samples. Because of the very small number of 
counts in individual decays (< 4 in r/6 and individual 
Xb.j transitions) statistically significant results for rji, and 
XbJ transitions can only be obtained by summing over 
all the decay channels to obtain ^iBiB^i- For the rjb, 
these are listed in Table I. Admittedly, these are still 
rather crude results, and statistical errors do not repre- 
sent the uncertainties reliably. A better measure of the 
uncertainty is provided by comparing the sum BiB2i 
over the three XbJ states and the 5 decays measured 
by us for each (-)• 3(7r+7r~), 4(7r+7r~), K+K~2{'k+-k~), 



') with the sum of the 



published results for them [8]. Our results, based on dif- 
ferent event selections and much smaller statistics, agree 
with the published results within a factor 1.5. 

Based on the fitted masses of the Xbj{^P) resonances, 



Sources of 


T(15) Results 


T(2S') Results 


oy OttJlllclllC UllCt;! LcLlllLico 


AMhf 




AAfhf 




(ranges of their variations) 


(MeV) 


(%) 


(MeV) 


(%) 


iNumDer oi i yTio j 




-1-9 




-1-9 


Mass Calibration (from Xb.i) 


±2.0 




±2.0 




Reconstruction & FID 




±16 




±16 


Detector Resolution (±10%) 


±0.2 


±1 


±0.1 


±5 


Fit Range 


±0.8 


±14 


±0.3 


±7 


12.5-20 250-300 MeV 










F(?76(15))=5-15 MeV 


±0.8 


±11 


±0.1 


±12 


r(776(25)) =2.5-7.5 MeV 










T{nS), 1st— 7-3rd order cxpon. 


±0.2 


±7 


±0.2 


±5 


Bin Size (1-2.5 MeV) 


±0.1 


±1 


±0.4 


±4 


Total 


±2.3 


±25 


±2.1 


±23 



we conservatively assign the systematic uncertainty of 
±2.0 MeV in our mass calibration. The systematic un- 
certainties in our results due to other possible sources 
were taken to be equal to the maximum variations in 
AMhf and B\ x B-i found in varying the source parame- 
ters from their nominal values. The range of variations 
for different parameters and the maximum variations ob- 
served in AMhf and B\ x S2 are listed in Table II. The 
uncertainty in number of T(nS') produced is estimated 
to be 2%. Event reconstruction and PID uncertainties 
vary for the different decay modes. The maximum un- 
certainties correspond to the largest multiplicity decays. 
To be conservative, in Table II we assign these maxi- 
mum uncertainties, 16% for r?h(15') and rib{2S), also to 
the sum of all modes. The systematic uncertainties added 
in quadrature are ±2.3 MeV in AMhf (15) and ±2.1 MeV 
in AMhf (25), and ±25% and ±23% in the corresponding 
Bi X B2- Our final results arc obtained by adding these 
in quadrature to the statistical errors in Table I. 

Our results for the IS state, M{rib{lS)) = 9393.2 ± 
4.1 MeV and AMhf (15)55 = 67.1 ± 4.1 MeV, agree with 
previous determinations [2]. Our identification of the 
r]b{2S) state at a ~ 5(7 level leads to 

M(77fc(25)) = 9974.6 ± 3.1 MeV 
AMhf (25)55 = 48.7 ± 3.1 MeV 

In summary, we have presented evidence for the first 
successful observation of r]b{2S), and the hyperfine split- 
ting of the bottomonium 2S state. Unqucnchcd lattice 
predictions for radial excitations are admittedly not yet 
very reliable [9-11]. The latest of these calculations [12] 
by the HPQCD CoUaboration obtains AMhf(15)55 = 
70 ± 9 MeV, in good agreement with its experimental 
value, and predicts AMhf (25)55 = 35 ± 3 MeV [12]. 

This investigation was done using CLEO data, and as 
members of the former CLEO Collaboration we thank it 
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ful comments. This research was supported by the U.S. 
Department of Energy. 
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